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Full counting statistics is a powerful tool to characterize the noise and correlations in transport 
through mesoscopic systems. In this work, we propose the concept of conditional spin counting 
statistics, i.e., the statistical fluctuations of a spin up (down) current given the observation of a 
given spin down (up) current. In the context of transport through a single quantum dot, it is 
demonstrated that a strong Coulomb interaction leads to a conditional spin counting statistics that 
exhibits a substantial change in comparison to that without Coulomb interaction. It thus can be 
served as an effective way to probe the Coulomb interactions in mesoscopic transport systems. 
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I. INTRODUCTION 

The exploration of full counting statistics (FCS) in 
mesoscopic systems has vital roles to play in gaining 
deep insight into microscopic mechanisms of transport 
and temporal correlations between charge carriers which 
are not available from conventional measurements of 
time-averaged current alone [H-Q- Particularly, recent 
advances in nanotechnology have made it possible to 
measure electron transport processes that take place at 
single-electron level p|-[l3|. All statistical cumulants of 
the number of transferred particles can now be extracted 
experimentally. 

Theoretical study of FCS based on the scattering ap- 
proach turns out to be very powerful for characterizing 
statistics of noninteracting electron transport through a 
variety of systems [3, [T3l - [l5| . Current fluctuations would 
satisfy the Poissonian process if the tunneling events 
were statistically independent. Normally, electric cur- 
rent noise is suppressed from the Poisson value due to the 
Pauli exclusion principle which forbids two electrons with 
the same quantum states to be superimposed [l^ [T7| . 
On the other hand, Coulomb interaction also serves as 
another important mechanism that can correlate wave 
packets. Particularly, with continued miniaturization of 
the system size, electron interaction effects become in- 
creasingly important in mesoscopic transport [isl [l9| . 
Recent study shows that the noise is indeed more sen- 
sitive to interactions than the average current [20l - l25| . 
Depending on different physical regimes concerned, the 
Coulomb repulsions may either decrease or increase noise 
correlations [26l - l32| . However, for a given mesoscopic sys- 
tem, the effects on the noise features due to Fermi statis- 
tics and that due to Coulomb interactions are normally 
intimately mixed. 

It is clear that the Pauli principle acts only on electrons 
with the same spin, while electrons with opposite spins 
are only correlated by the Coulomb interaction. It thus 



inspires us to propose in this work the concept of con- 
ditional spin counting statistics: The statistical fluctua- 
tions of a spin up (down) current given the observation 
of a given spin down (up) current. The theory is applied 
to a simple system-a single quantum dot (QD) tunnel- 
coupled to two normal electrodes. Although the net cur- 
rent is not spin polarized, the transport of up and down 
spins are actually correlated to each other via Coulomb 
interaction. It is demonstrated that the strong Coulomb 
interaction leads to a conditional spin counting statis- 
tics that undergoes a substantial change in comparison 
to that without Coulomb interaction. It thus may offer 
an effective and transparent way to probe the Coulomb 
interactions in various mesoscopic transport systems. 

The paper is organized as follows. In section [Hi we 
describe the single QD system in different bias configu- 
rations, such that the effectiveness of the Coulomb in- 
teractions can be fully taken into account. Section IIIII 
is dedicated to the theory of conditional spin counting 
statistics. Application of the theory to QD system is 
demonstrated in Section IIVI with focus on its effective- 
ness in probing Coulomb interaction. It is then followed 
by the conclusion in section FVl 



II. THE MODEL 

We consider electron transport through a single QD 
with Coulomb interaction, as schematically shown in 
Fig.H The entire system is described by the Hamilto- 
nian H = Hb + Hqu + H' , with 



a— L,R k(T 
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Here H-q models the noninteracting electrons in the left 
(a=L) and right (a=R) electrodes, with c^j,^ {caka) the 
electron creation (annihilation) operator in the corre- 
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(a) Model (I) 




(b) Model (II) 
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FIG. 1: Schematic setup of transport through a single quan- 
tum dot in different bias configurations. Model (I): The bias 
is large enough to overcome the Coulomb interaction. Model 
(II): A small bias is applied across the QD. The Coulomb re- 
pulsion allows at most one electron (either spin up or spin 
down) to be occupied. In both cases, the Fermi functions can 
be approximated by either one or zero, so the temperature is 
not involved here. 



sponding electrode. The electron distributions in the 
electrodes are governed by the electrochemical potentials 
/iL and /iR. In this simple model, we choose a symmet- 
ric bias voltage, such that — eV/2 and ji-R. = —eV/2. 
Hqd describes the QD with one spin-degenerate energy 
level and the Coulomb interaction U on the dot, where 
Tia = d^dcr is the occupation operator, with dj. (d^) the 
electron annihilation (creation) operator in the QD. Elec- 
tron tunneling between electrodes and QD is depicted by 
H' . The tunnel-coupling strength is characterized by the 
intrinsic tunneling width Ta{Lij) = 27r(7Q(a;)|tQ[^, where 
ga(y>) is the density of states in the electrode a. Here- 
after we assume flat bands in the electrodes, which yields 
energy-independent couplings Fq,. Throughout this work, 
we set fi, = 1 for the Planck constant, unless stated oth- 
erwise. 

By specifying which excitation energies lie within the 
energy window defined by the Fermi levels /zl and /zr, the 
following two configurations will be considered. Model 

(i) : The bias is large enough to overcome the Coulomb 
interaction, as schematically shown in Fig.[TJa). Model 

(ii) : Only the single level £-|-/j^ is within the bias, as shown 
in Fig.[TJb). The charge transport is strongly correlated. 
The present model is simple; however, we will show that 
it is adequate to address the essence of the conditional 
spin counting statistics and its effectiveness in probing 
Coulomb interaction. 



III. FULL COUNTING STATISTICS VERSUS 
CONDITIONAL SPIN COUNTING STATISTICS 

The FCS of charge transfer can be described in the 
framework of counting field-dressed approach, in which 
electron jumps between reduced system and electrodes 
are taken into account by introducing a corresponding 



counting field x- This results in a x-resolved quantum 
master equation, which in the Born-Markov limit can be 



formally written as [SSHSTf ^p(x,i) = C^p{x,t). The 
cumulant generating function (CGF) is determined from 
the lowest eigenvalue of C^^. For instance, the CGF for 
model (I) is given by [s^l 



Hx) = -tc{TL + Tr - v/(rL-rR)2+4FLFRe-ix}, 

(2a) 

while for model (II), the CGF reads 



= -^{2rL + Tr - v/(2FL-FR)2+8rLrRe-ix}, 

(2b) 

where tc is the counting time. Apart from an effec- 
tive doubling of the rate Fl, the result for model (II) 
[Eq. ()2bp ] is qualitatively similar to that of model (I) 
[Eq. ([2a| ]. Therefore the charge counting statistics is not 
the best possible probe of interactions. We will show 
later, on the contrary, the behavior of conditional spin 
current noise is completely different, and can be served 
as a sensitive diagnostic tool to probe the Coulomb in- 
teraction. 

Now we introduce the concept of conditional spin 
counting statistics: The statistical fluctuations of the 
spin-t (-1) current given the observation of a given spin- J, 
(t) current. These statistics may be calculated from the 
conditional distribution functions, P(/^|/^) or P(/^|/-|-), 
the probability of observing one spin current component 
given an observation of the other. To obtain the mixed 
generating functions of conditional spin counting statis- 
tics, we introduce spin-resolved counting fields ^^id 
Xi, which are used to characterize, respectively, jumps of 
up and down spins through a specific junction. The cor- 
responding (x-|-,x^)-resolved master equation can be for- 
mally written as ■§^p{x^,Xl,t) = M(xt, Xi)p(xt' X;, 
In the stationary limit, the counting time tc is much 
longer than the time of tunneling through the system, the 
joint generating function -F(xt7 Xi) ^-nd consequently the 
counting statistics is determined from the minimal eigen- 
value of M according to 



J^iXtiXl) = -Amin(XT'X;)ic 



(3) 



where Amin satisfies A 



min|xt^-0,X4-->0 



0. 



From the definition of the joint generating function, the 
joint probability distribution of transmitted spins can be 
extracted by Fourier transforming on both variables. 



P(7Vt,7V;,tc) 



(2^)2 



(4) 



Replacing A'^ ~ I^tc and — I^tc gives the joint prob- 
ability distribution of the two currents. In the stationary 
limit (tc 3> rL^,rp"^), it is justified to evaluate the in- 
tegral in the saddle point approximation (37l - l40| . The 
dominant term contributing to the joint distribution is 
then given by 



{Amin + i/tXT + i^iX;}- (5) 
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The mixed generating functions of conditional spin 
counting statistics may be calculated by only Fourier 
transforming on one of the above variables in Eq. ([5]) . 
For instance, integrating over the xt variable yields 
■^ixt'-^l)- to the Bayes's theorem, which relates 

the joint distribution and conditional distribution func- 
tions P(/t,/;) - P(/t|/;)P(/;)^ = P(/;|/t)P(/t), the 
mixed generating function then is given by J-(x-t\Ii) — 
J-{xti ^i) ~ -^(0' ^i)- The statistical fluctuations of spin 
t current given the observation of a given spin J, current 
can be obtained from T{x'\-\Ii) by performing derivatives 
with respect to the counting field 



(It. 



(-ia 



/ L - 



x-r) 



--^(Xtl^i)lxt^o 



(6) 
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Analogously, the conditional spin counting statistics of 
spin ], current can be evaluated from J^(x^|/-t-)- 



FIG. 2: Conditional spin t current and its noise versus spin 
4, current. The dashed and solid lines correspond to results 
of the model (I) and (II), respectively. We use Fl = F as the 
energy unit, and set rR=4r. 



IV. RESULTS AND DISCUSSION 



First let us focus on the model (I), as schematically 
shown in Fig.[IJa). The involving states are |0)-cmpty 
dot, I t/i)-occupation by a spin-t/J, electron, and |d)- 
doubly occupied. The X4^)-resolved quantum master 

equation reads ■§^p{xt,Xl) = Mi{xt,Xl)p{xt,Xi), where 
Ml is a 4x4 matrix 



Ml 



( -2Fl 
\ 



-Fl - Fr 






-Fl - Fr 



.0 \ 

eiX4,rR 

e'xtrR 

-2Fr ; 



(7) 



The minimal eigenvalue of M\ is then determined, which 
leads to the joint distribution Pi(/-|-, /j^, tc) of spin t and 
J, currents 



log Pi 



<T=t4 



, (8) 



with = 2/^ + v/4/2 + (Fl - Fr)2. Apparently, 
the joint probability factorizes Pi{Ii-,Ii) = P{Ii-)P{Ii), 
which implies that spin t and ]. currents arc uncorrelated. 
The resultant mixed generating function for conditional 
spin counting statistics then reads 



^lixtlh) 



tr 



{ri 



'(rL-FR)2+4FLFRe-ixt}. 

(9) 

which is the same as the charge counting statistics for 
model (I) in Eq. ([2a| , except for an overall factor of | . 
This is because the spin transport through the QD in- 
dependently, and each spin component contribution 50% 
of the total current. The cumulants of conditional spin 
counting statistics are thus the same as those of uncon- 
ditional one. For instance, the first and second mixed 



cumulants are given, respectively, by 

Tl+Tr' 



2e(/t)/, 2(Fl 



(10a) 
(10b) 



In this case the system can be mapped onto that trans- 
port through a single level without Coulomb interaction. 
One may apply analogous analysis to the conditional 
counting statistical of the spin J, current given the ob- 
servation of a spin ^ current, and the same results as 
Eq. will be obtained. 

Fet us now consider the conditional spin counting 
statistics for model (II) as shown in Fig.[TJb), where at 
most one electron can be occupied and the spin currents 
are strongly correlated. Thus the involving states are re- 
duced to |0)-empty dot, |t/4.)-occupied by a spin t/4- elec- 
tron. The spin-resolved quantum master equation reads 
^p(xt'Xi) = ^ii(xt:X;)p(xt,Xi): where the 3x3 rate 
matrix is given by 



-2Fl e'^^TTR e'x^Fr 
Mtt = I Fl -Fr 
Tl -Fr 



(11) 



The minimal eigenvalue can be readily evaluated, which 
results in the joint probability for the spin currents for 
the model (II) 



logPii 

tc. 



=i(A-2FL-FR)-^/,log 



A/^ 
FlFr 



(12) 



where A = 2(7^ + 4) + v/4(/t + /i)^ + (2Fl - Fr)^. The 
term inside the logarithm shows unambiguously that the 
spin t and currents arc correlated. Further integrating 
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FIG. 3: Schematic setup for measuring the conditional spin 
counting statistics by employing a device made up of four 
spin-polarized terminals (see text). 



over the xt variable gives rise to the mixed generating 
function 



-^ii(xtl4) = | { C?(xt)-G(0)-2/^log 



2/^ + G(xt) 



24 + G(0) 



(13) 
The 



with Gixt) = y^4/2 + 4e>^trLrR + (2rL-rR,) 

conditional spin counting statistics may be calculated via 
taking derivatives with respective to x^- For instance, the 
first mixed cumulant yields the conditional spin current 



{It)i 



FtF 



Li R 



2h + GiO)-' 



(14a) 



and the second mixed cumulant corresponds to the con- 
ditional shot noise of spin t current 



2e(/t)i 



2rLrR \ 

[2/^ + G(0)]G(0)y'' 



(14b) 



to zero for ^ 00. It might be larger or smaller than 
that of model (I), depending on tunneling rates Fl/r 



and I^. For a given spin J, current 



the 



. . 2rL+rR ■ 

conditional spin current and the unconditional one co- 
incide ((/-f)/^ = (-^t) uncond), which corrcspouds to the 
most likely spin current configuration for non-polarized 
electrodes. The conditional spin current noise for the 
model (II) however increases with rising I^. It reaches 



the maximum i 



as I\ 



00, which reflects the rare tun- 



neling events of up spins through the QD system given 
an extremely large spin ], current. These unique noise be- 
havior together with their striking distinction compared 
to the model (I) thus can be served as a sensitive tool to 
probe the Coulomb interaction in mesoscopic transport. 

Finally, let us propose a setup for the measurement 
of conditional spin counting statistics by considering a 
four-terminal device j4ll |43 |. as shown in Fig.|31 The 
four electrodes are made up of fully spin polarized ferro- 
magnetic metals. The left two electrodes (LI and L2) of 
opposite spin polarizations arc kept at the same chem- 
ical potential //l- Likewise, electrodes Rl and R2 with 
opposite polarizations on the right have the same chem- 
ical potential ^r . If the junction parameters for LI and 
L2 are the same, and those for Rl and R2 are identical 
as well, the net current transport through the QD is not 
spin polarized. This setup thus can be mapped onto the 
model that we have analyzed. In spite of non-polarized 
current, it is possible to measure separately the up or 
down spin current in each of the four electrodes, which 
eventually enables access to the conditional spin counting 
statistics. We expect that these predictions can be tested 
in quantum transport experiments in the near future. 



where G(0) = ^Alf + ATI + ^r- Both show a radical 
change in comparison to those without Coulomb interac- 
tion, sec Eq. (|10|) for the model (I). One may also compare 
the conditional spin current and noise with the uncondi- 
tional quantities [41 1 



FtT 



Ll R 



2e(/t) 



uncond 



uncond 

1 

" 2 



2Fl + Fr' 
2FLrR 



(2rL 



rR)2 



(15a) 
(15b) 



which are similar to the conditional cumulants without 
Coulomb interaction [Eq. (fTO|) ] , except for an effective 
doubling of Fl. Thus, the conditional spin counting 
statistics manifests itself unambiguously as a sensitive 
probe of Coulomb interaction. 

The numerical results for the conditional spin t current 
and its noise versus spin 4, current are plotted in Fig.[2Ka) 
and (b) , respectively. For the model (I) , the spin '[ and 4, 
currents are uncorrelatcd, thus the conditional spin f cur- 
rent and its noise [Eq. ([T0|) ] are independent of the spin 
I current, as shown by the dashed lines in Fig. [21 The 
conditional spin t current (/^)/^ for the model (II) de- 
creases monotonically with the spin J, current, and tends 



V. CONCLUSION 

In summary, we have proposed the concept of con- 
ditional spin counting statistics: The statistical fluctu- 
ations of a spin up (down) current given the observa- 
tion of a given spin down (up) current. By applying 
the theory to a simple system-a single quantum dot, we 
demonstrated that in the presence of a strong Coulomb 
interaction the conditional spin counting statistics under- 
goes a substantial change in comparison to that without 
Coulomb interaction. It thus can be served as an effec- 
tive way to probe the Coulomb interactions in mesoscopic 
transport systems. Furthermore, the present theory may 
also open the possibilities to get a deep understanding of 
the spin-resolved bunching behavior and its connection 
to the noise characteristics [ij. 
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